Pseudomodes as an effective description of memory 
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We investigate the non - Markovian dynamics of two -state systems in structured reservoirs. We 
establish a connection between two theoretical quantum approaches, the pseudomodes [B. M. Gar- 
raway, Phys. Rev. A 55, 2290 (1997)] and the recently developed non - Markovian quantum jump 
method [J. Piilo et al, Phys. Rev. Lett. 100, 180402 (2008)]. This connection provides a clear 
physical picture of how the structured reservoir affects the system dynamics, suggesting the role of 
the pseudomodes as effective description of environmental memory. 
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I. INTRODUCTION 

Recent advances in experimental techniques for coher- 
ent control of small quantum systems have paved the way 
to a series of spectacular experiments aimed at both test- 
ing fundamental features of quantum theory and imple- 
menting logic gates for quantum information processing 
[TJ |2]. Quantum properties, however, are very fragile. 
Any interaction between quantum systems and their sur- 
roundings gives rise to decoherence and dissipation phe- 
nomena, destroying the quantumness of the state of the 
system. For this reason during the last decade several 
theoretical and experimental studies have been devoted 
to the investigation of the dynamics of open quantum 
systems [3j El |5] . 

One of the approaches to the description of open quan- 
tum systems consists in separating the total system into 
two parts, the quantum system of interest and the sur- 
rounding environment. The environment is often mod- 
eled as an infinite collection of quantum harmonic oscil- 
lators in thermal equilibrium [3|. This type of reservoir 
describes, e.g., the quantized electromagnetic field. One 
of the key quantities characterizing the reservoir, and 
therefore determining the open system dynamics, is the 
spectral distribution or structure function. This quantity 
describes the frequency - dependent coupling between the 
system and the continuum of harmonic oscillators form- 
ing the environment. 

In many physical situations the system - reservoir cou- 
pling strength does not depend strongly on the frequency 
of the reservoir oscillators. In this case the reservoir 
spectral density can be conveniently approximated by a 
flat spectrum. One typically refers to such systems as 
Markovian open quantum systems. A typical feature of 
the dynamics of Markovian open quantum system is the 
irreversible flow of energy and/or information from the 
system to the environment. 

In certain physical contexts, however, the quantum 
system of interest interacts with "structured" reservoirs, 
whose spectral density strongly varies with frequency. 
We refer to such systems as non - Markovian open quan- 
tum systems. Non - Markovian dynamics is characterized 
by the existence of a memory time scale during which 



some energy/information that has been transferred from 
the system to the environment feeds back into the sys- 
tem. 

Many solid - state systems such as Josephson junctions, 
display strongly non - Markovian dynamics [5 . Moreover, 
often the e.m. field surrounding a quantum system, e.g. 
an atom, can be conveniently engineered to prevent or 
inhibit the occurrence of decoherent processes such as 
atomic spontaneous emission. Atoms placed inside cavi- 
ties or photonic band gap materials are examples of quan- 
tum systems interacting with such engineered structured 
reservoirs [B] . The formalism we develop here for a two - 
level atom can be applied to any two - level system in 
bosonic structured reservoir, e.g, NV vacancy in diamond 
embedded in photonic band gap [TJ. 

The theoretical description of the dynamics of non- 
Markovian quantum system is usually very complicated, 
only a few simple systems are amenable to an exact so- 
lution [3] . A number of methods have been formulated 
for treating non-Markovian dynamics [8j [9j [10] but the 
connection between these methods has so far remained 
unexplored. Most importantly, due to the mathemati- 
cal difficulties in dealing with non-Markovian systems, 
a simple intuitive physical picture of the memory of a 
non - Markovian reservoir and of how such memory al- 
lows to partly restore some of the coherence lost to the 
environment is highly desirable. Such a simple effec- 
tive description is our main result. By connecting two 
non-Markovian approaches, the pseudomodes [8] and the 
non-Markovian quantum jumps (NMQJ) methods [5], 
we are able to identify, for simple exemplary cases, where 
the memory of a non-Markovian reservoir resides. 

Our results provide insight into the problem of the exis- 
tence of pure state trajectories and the physical meaning 
of the master equation unravelling for non - Markovian 
systems. This issue has been recently considered in Rcfs. 
[TT1 112] . Both the pseudomodes and the NMQJ ap- 
proaches are consistent with the interpretation given in 
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II. PSEUDOMODE AND NMQJ METHODS 

A. Time-local master equation 

Consider a two - level atom interacting with a struc- 
tured electromagnetic reservoir in the vacuum state. The 
Hamiltonian of such a system in the rotating wave ap- 
proximation is 



H = h0J()(7j 



-2>.-t 



xa{ax + ^2(hglax<T + +R.c.) (1) 



where a± are the Pauli raising and lowering operators 
for the two - level system, luq is the atomic transition fre- 
quency, a\ and a x are the annihilation and creation op- 
erators for the mode A of the field having frequency ui\ 
and coupling constant g\. We assume that initially only 
one excitation is present in the system. The state of the 
total system at time t takes the form 

\m) = colff, A ) + d{t)\e, Ox) + c A(*)b, ( 2 ) 



with \g, Ox), |e, Ox), and \g, l\) the states containing zero 
excitations, one atomic excitation and one excitation in 
the A-mode of the e.m. field, respectively. 

The problem of a two -level system interacting with 
a zero temperature reservoir is in principle exactly solv- 
able using Laplace transforms. The exact non-Markovian 
master equation describing the dynamics of the atomic 
system takes the form [3J 

~~T7~ = ^-[v+v-,Pa\ +j(t)[a-PA<r+ - \{<r+a-,p A }], 
at li 2 

(3) 

where pa is the atomic density operator. The time de- 
pendent Lamb -shift S(t) and the time dependent decay 
rate j(t) are given by 

The master equation ^ can be simulated by means 
of the non-Markovian quantum jump (NMQJ) method 
[9] which extends the Monte Carlo wave function ap- 
proach JT3] to non - Markovian systems with negative de- 
cay rates. 



B. Pseudomode method 

Alternatively, one can investigate the dynamics using 
the pseudomode theory [H [14]. This method relies on 
the strong connection between the atom dynamics and 
the shape of the reservoir spectral distribution. More pre- 
cisely, the key quantities influencing the time evolution of 
the atom are the poles of the spectral distribution in the 
lower half complex ljx plane. By introducing some aux- 
iliary variables, called pseudomodes, defined in terms of 



the position and of the residue of the poles of the spectral 
distribution, one can derive a Markovian master equation 
in the Lindblad form for the extended system compris- 
ing the atom and the pseudomodes. This exact master 
equation describes the coherent interaction between the 
atom and the pseudomodes in presence of decay of the 
pseudomodes due to the interaction with a Markovian 
reservoir. 

For a Lorentzian spectral distribution the pseudomode 
approach leads to the following master equation 



dp 
~dl 



Y, 



[Hq, p] — — [a}ap — 2apd* + pa' a 



(5) 



where 



Hq = L}Qa + <j- + La c a) a + Q,o[a '<r_ + a<7+] (6) 



and p is the density operator for the atom and the pseu- 
domode. Since a Lorentzian function has only one pole 
in the lower half complex plane, the atom interacts with 
one pseudomode only as displayed in Fig. 0(a). The con- 
stants uj c and r are respectively the oscillation frequency 
and the decay rate of the pseudomode and they depend 
on the position of the pole Z\ = to c — iT/2 while Slo is the 
pseudomode coupling constant. 

Both master equations Q and §5§ are exact. Hence 
we expect to obtain an equation of motion for the atom 
of the form of Eq. (|3| by tracing out the pseudomode in 
Eq. ([5]). Additionally it is interesting to see the expres- 
sions of the coefficients 7(f) and S(t) as functions of the 
pseudomode amplitude. Indeed the equation obtained 
from Eq. ^ has the form 

dp a Ait) 1 

-jr = -of [v+v-iPAl +B{t)[(T-pA(T+ - -{a + <j-,p A \\, 

(7) 



where 



and 



A(t) = 2 



uj + Q 



Re{ Cl (t)bl(t)} 
IciWI 2 



B{t) = 2ft, 



Im{ Cl (t)bl(t)} 
' IciWI 2 ' 



(8) 



(9) 



where bi(t) is the pseudomode amplitude. By using the 
differential equations governing the atom - pseudomode 
dynamics [5] 

i— ci = cj ci + floh, i—bi = zibi+fl ci, (10) 
it is easy to prove that A(t) = S(t) and B{t) = j(t). 



C. The connection between pseudomode and 
NMQJ approaches 

Once we have proven the equivalence between the 
two master equations, we focus on the simple case of 
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a Lorentzian spectral distribution off- resonant with the 
atomic transition frequency (damped Jaynes- Cummings 
model with detuning). 

In the strong coupling regime T <C the atomic dy- 
namics is strongly non - Markovian. The excited state 
population of the atom oscillates in time indicating that 
the energy dissipated into the environment flows back 
into the system as a consequence of the reservoir mem- 
ory. At the same time the atomic decay rate 7(f) attains 
negative values. The key role of the pseudomode is ex- 
posed when we look at the time derivative of the pseu- 
domode population and, using the differential equations 
in Eq. (10 1, we obtain 



dt 



•r|6i(i)| 2 = 7(*)|ci(i)| S 



(11) 



The equation above shows that the compensated rate of 
change of the pseudomode population, given by the left 
hand side of Eq. (Ill (where the effect of the pseudomode 
leakage is removed), is directly related to the atomic de- 
cay rate. This means that the information about the 
dissipative dynamics of the atom into the structured 
reservoir are all contained in the pseudomode dynam- 



ics. Equation (111 and its physical interpretation is one 



of the main results of the paper. 

We conclude this section analyzing the connection be- 
tween the NMQJ unraveling of the master equation in 
Eq. @ and the Monte-Carlo wave function unraveling of 
the pseudomode master equation in Eqs. (|5| and (21 1. 
In the NMQJ description the ensemble members living 
in the Hilbert space of the system are always in a pure 
state. In particular the density matrix of the ensemble 
can be written as 



PAit) 



Np(t) 
N 



\Mt)HMt)\ 



Ni(t) 
N 



IV-iXM (12) 



where 



\Mt)) = C g {t)\g) + C e {t)\e) 
Hi) = Iff), 



(13) 
(14) 



N is the total number of ensemble members, Ni(t) is the 
number of ensemble members who have jumped into the 
ground state, and No(t) is the number of members who 
have not jumped, or which have gone through a jump- 
reverse-jump cycle due to the negative decay rate [S]. 

In the pseudomode description the unraveling of the 
master equation (|5| is in the extended Hilbert space con- 
taining the pseudomode. The density matrix of such en- 
semble, expressed in the atom-pseudomode basis, is the 
following: 



N p 



l h 



(15) 



where 



\<(t)) = C^(t)\g,0)+C^(t)\9A)+C^ (t)\e,0)(16) 

\tf) = \g,o), (17) 



N p is the total number of ensemble members, N p (t) is 
the number of ensemble members who have decayed into 
the atom-pseudomode ground state via a pseudomode 
jump, and N p (t) is the number of members who have 
not jumped. If we want to look at the time evolution of 
the ensemble members in the atomic Hilbert space only, 
we have to trace out the pseudomode auxiliary degree 
of freedom. This leads to the following reduced atomic 
density matrix: 



PA(t) 



■■^((K (t)\ 2 + \Cm 2 )\9){9\ 



+ |Cf (i)| 2 |e)( e 
+C^{t)C^{t)\g){e 



C^(t)C^(t)\e){g\ 

N 1 P (t) l 
N p 



(18) 



\g){g\ 



in which the ensemble members are clearly in a mixed 
state. A comparison between Eqs. (12 1 and (18) illus- 
trates the connection between the two unravelings. In 
particular it is illustrative to consider the ground state 
population, 



(g\ P A(t)\g) = 



Ni(t) , N (t) 



N 
N p 



N 
N P (t) 
N p 



c 9 (t)\ 2 



(\C P (t)\< 



|c£(t)| 2 ), 



(19) 



further showing the unravelings connection. 



D. The pseudomode as effective memory 



The interpretation of Eq. (11 1 is particularly interest 



ing in the light of the NMQJ method, where negative de- 
cay rates lead to reversed quantum jumps. Consider, for 
example, an atom initially prepared in a generic superpo- 
sition of ground and excited state performing a quantum 
jump to its ground state at a certain time t' . If the decay 
rate *y(t) becomes negative at t > t' , the superpositions 
destroyed by the earlier normal jump can be restored by 
a reversed jump. In fact a reversed jump takes the atom 
to the state into which it would have evolved if the previ- 
ous quantum jump had not occurred. Thus in the NMQJ 
framework one can characterize the period of negativity 
of the decay rate as the period of time in which memory 
effects and restoration of quantum superpositions occur 
through reversed quantum jumps. Stated another way, 
the reverse jumps describe the process through which the 
system recovers part of the information that leaked into 
the environment. This is confirmed by the dynamics of 
the von Neumann entropy of the atom which shows an os- 
cillatory behavior following the oscillations of the atomic 
decay rate "f(t). In particular the atomic von Neumann 
entropy decreases when the decay rate is negative indi- 
cating a temporary reduction of the mixedness of the 
atomic state and a recovery of coherence. 

Equation (111 states that if the decay rate is negative 



then the compensated rate of change of the pseudomode 
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FIG. 2: The solid line is the decay rate for a two -level atom 
in a Lorentzian structured reservoir. The dashed line is the 
compensated rate of change of the pseudomode population 
d\bi{t)\ 2 /dt + r\bi(t)\ 2 . The units for the rates on the ver- 
tical and horizontal axis are 70 and I/70 with 70 = 4f2o/r 
Markovian decay rate of the atom. We have taken the values 
T = 0.6 70, fio = \/0.15 70 and lo c — loq = 4T. 



FIG. 1: Diagramatic representation of the atom -pseudomode 
dynamics, (a) Atom interacting with a Lorentzian structured 
reservoir: the atom interacts with a single pseudomode which 
leaks into a Markovian reservoir, (b) Atom in a "simple" pho- 
tonic band gap: we see a more complex memory architecture, 
the second pseudomode acts as a memory for the atom, while 
the first pseudomode acts as a memory for the first one. 



population is negative as well, as clearly shown by Fig. [2] 
So whenever the atom increases its excited state popula- 
tion the pseudomode must deplete. This equation, there- 
fore, establishes a link between the restoration of coher- 
ence, typical of a reversed jump, and the pseudomode 
depletion. This observation suggests an interpretation of 
the pseudomode as that part of the reservoir from which 
the atom receives back information and probability due 
to memory effects. In other words the pseudomode can 
be seen as an effective description of the reservoir mem- 
ory. This is further shown by the dynamics of the mutual 
information f2\ [TS] between atom and pseudomode which 
perfectly follows the oscillation of the pseudomode popu- 
lation. However the pseudomode is not a perfect storage 
place since the efficiency of the information restoration 
from the pseudomode to the atom depends on the pseu- 
domode loss rate T. 



III. MEMORY STRUCTURE AND 
PSEUDOMODE ARCHITECTURE 



We now generalize our study to a more complicated 
spectral distribution, namely the inverted Lorentzian mo- 



del describing in a simplified way photonic band gaps 



D(u) 



11 1 



W 2 T 



2-L 2 



(w-wc) 2 + Cri/2) 2 (^-^ c ) 2 + (r 2 /2)2- 

(20) 

The negative Lorentzian introduces a dip into the den- 
sity of states leading to the inhibition of atomic sponta- 
neous emission in the region of the dip. In particular, for 
Wx/Tt = W2/T2, the spectral distribution in Eq. ((20 1 
presents a perfect gap, D{uj c ) = 0. The exact pseudo- 
mode master equation is given by [8] 



dp 



where 



-%[H ,p] 

r" 

^-[ a l a 2P - 2a 2 /9<4 + pa\a 2 ], 



2a\pa\ + pa\ai] 



(21) 



Hq =ujo<t + (T- + uj c a\ai + LU c a\a2 + ilo[c4' J - 
+ V(a\a 2 + ai(4), 



CJ2C4 



(22) 



where a\ and a 2 are the annihilation operators of the two 
pseudomodes decaying with decay rates r'j = I4 y ir 2 — 
W^ri and r' 2 = W x Tx - W 2 T 2 , respectively. The two 
pseudomodes are coupled and V = t/WiW2(Ti — T 2 )/2 
is the strength of the coupling. Figure [l] (b) shows the 
atom - pseudomodes architecture in this case. The atom 
interacts coherently with the second pseudomode, which 
is in turn coupled to the first one. Both pseudomodes are 
leaking into independent Markovian reservoirs. In the 
case of a perfect gap only the second pseudomode leaks. 
The set of ordinary differential equations associated to 
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the master equation (21) is 



IV. CONCLUSION 



*37 c i = ^oci + ^002, 
at 

■ d 

i— ai = z x a\ + Va 2 , 

i^a 2 = z' 2 a 2 + Vai + QqCi, 



(23) 



where c\, a\ and a 2 are the complex amplitudes for the 
states with one excitation in the atom, one excitation in 
the first pseudomode, and one excitation in the second 
pseudomode, respectively. The position of the true poles 



is z\ 



and z' 2 = uj c - iT' 2 /2. 



Similarly to the calculations for the Lorentzian case 
one can show that, after tracing out the two pseudomodes 
in Eq. (21), and with the help of Eqs. (23 1, one obtains 
the non - Markovian master equation (|7j) for the atom, 
where A(t) and B(t) are given by the same expressions 
in Eq. |8]) and |9]) provided that we replace b\{t) with 
oa(t). 

The non - Markovian dynamics of the atom is linked 
to the coherent variation of both pseudomodes by the 
following equation 



d\ ai {t)\^ 



dt 



+rik(i)| 2 



d\a 2 (t)\< 



dt 



-n|a 2 (t)| 2 = 7 (i)| Cl (t)| 2 . 

_ (24) 
This equation generalizes Eq. (Ill to the more com- 



plex reservoir structure considered here. Moreover, using 
Eqs. (23 1 we obtain a relation connecting the dynamics 



of the first and second pseudomodes, 



d\ai(t)\< 
dt 



-riMt)| s 



Jm{a 2 {t)a 1 (t)} 2 

lV 1 — 77VT2 KM*) • 

\a 2 (t)\ 2 



(25) 



Having in mind that -y(t) = 2fl lm{c 1 (t)b* 1 (t)} /\ci(t) | 2 
one sees that Eq. ( 25 1 has the same structure as Eq. (ITT| . 



In fact the compensated rate of change of the first pseu- 
domode population equals a time dependent coefficient 
times the second pseudomode population. Therefore the 
first pseudomode acts as a memory storage for the sec- 
ond one. Such a storage of memory of the second pseu- 
domode into the first one is perfect in the case of perfect 
gap where = 0. 



Both the pseudomodes and the reversed quantum 
jumps describe the non - Markovian backaction of a reser- 
voir on a quantum system. Our results show that the pic- 
ture of open quantum systems as comprised of a system 
and environment is a commonly held view that does not 
reflect the dynamics when the coupling is strong; the en- 
vironment divides into memory and non-memory parts. 
Establishing a connection in the case of a general struc- 
tured reservoir, i.e., for a generic form of the spectrum is 
an extremely complicated issue because both the number 
of pseudomodes and their equations of motion depend on 
the details of the reservoir structure. In this sense our 
results constitutes a first step in the direction of an un- 
derstanding of memory in non-Markovian systems. 

It is likely that a measurement on the reservoir will 
in general affect any memory that it contains, and so 
far it has not been possible to attach any measurement 
scheme to non-Markovian reservoirs, as Ref. [H] demon- 
strates. However, since we have seen that it is possible 
to make a separation of the reservoir into a memory part 
and a Markovian reservoir part, there remains the pos- 
sibility that we could make a measurement scheme that 
would address only the Markovian part of the reservoir. 
In the cavity-atom situation, this might correspond to a 
measurement on outside-modes, which have a significant 
amplitude outside the cavity. In such a case it is clear 
that the reservoir-memory for the atom actually resides 
in inaccessible inside mode, or pseudomode, belonging to 
the cavity. 
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